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ABSTRACT 

In this article, we explore the possibility of using the properties of gamma ray bursts (GRBs) 
to probe the physical conditions in the epochs prior to reionization. The redshift distribution of 
GRBs is modelled using the Press-Schechter formalism with an assumption that they follow 
the cosmic star formation history. We reproduce the observed star formation rate obtained 
from galaxies in the redshift range < z < 5, as well as the redshift distribution of the GRBs 
inferred from the luminosity-variability correlation of the burst light curve. We show that the 
fraction of GRBs at high redshifts, whose afterglows cannot be observed in R and I band due 
to Hi Gunn Peterson optical depth can, at the most, account for one third of the dark GRBs. 
The observed redshift distribution of GRBs, with much less scatter than the one available 
today, can put stringent constraints on the epoch of reionization and the nature of gas cooling 
in the epochs prior to reionization. 
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1 INTRODUCTION 

The usefulness of gamma ray bursts (GRBs) to probe the evolution 
of the universe is realized ever since the redshift measurements of 



& Ostriker (1997), Chiu & Ostriker (2000)). These studies have 



GRBs became possible (Kulkarni et al. 19981. A direct connection 
between the GRBs and star formation rate (SFR) in the host galaxy 
is expected in the two popular models for GRB progenitors (col- 
lapsar and binary coalescence). There are also growing observa- 
tional indications for t he association of sup e rnova-like compone nts 
in several afterglows (Galama et al. 1999; Bloom et al. 1999j ). If 
indeed the rate of formation of GRBs is related to the SFR of host 
galaxy, they can be use d as an effective probe o f the star formation 
histo ry of the universe ([Lamb & Reichart 2000|; Porcia ni & Madau 
2001; piardi & Loeb 200$ |Blain & Natarajan 2000} Bromm & 
Loeb 2002). Thus the redshift distribution of GRBs can comple- 
ment the studies based on Lyman break galaxies in understanding 
the cosmic star formation history. At present only few tens of GRBs 
have redshift measurement and the list is, by no means, statistically 
complete. Recently it has been suggested that the redshift distribu- 
tion of the GRBs can be probed using the global correlation that 
seems to exist between different observables of the afterglow. The 
redshift distribution of GRBs, albeit with large scatter, can be ob- 
tain ed using the correlation between GRB luminosity and v ariabil- 
ity ( Fenimore & Ramirez-Ruiz 2000 ; Reichart et al. 200U L loyd- 
Ronning, Fryer, & Ramirez-Ruiz 2002). 

The formation of galaxies and cosmic star formation histories 
have been studied extensively through a nalytical calculations an d 
numerical simulations (see, for example, Rees & Ostriker (1977 ), 
White & Rees (1978;), |Katz, Weinberg, & Hernquist (1996|), Gnedin 



been extended to investigate the cosmological implications of the 
GRBs. In this work, we take a simple model of the cosmic star 
formation history, based on Press-Schechter theory of collapsing 
haloes (in this work, the term 'halo' is used for referring to the dark 
matter halo) and elementary ideas on cooling and galaxy formation. 
The SFR density at low redshift in the model i s constrained by the 



rest-UV luminosity observations of galaxies (Treyer et al. 1998 
Cowie, Songaila, & Barger 1999j |Steidel et al. 1999|; S omerville, 



Primack, & Faber 2001). However, at high redshifts (say z > 6) 
the nature of reionization (viz. the epoch of reionization and nature 
of baryonic cooling in haloes) strongly influence the star forma- 
tion history. The main aim of this article is to use observed redshift 
distribution of GRBs to probe the physical processes in the pre- 
reionization epoch. Indeed, one cannot probe this epoch directly 
in the optical wave band due to high Gunn-Peterson optical depth 
even if the intergalactic medium (IGM) is not completely neutral. 



2 ANALYTICAL FORMALISM AND MODEL 
PARAMETERS 

We use the Press-Schechter (PS) formalism to obtain the co- 
moving number density of collapsed objects having mass in the 
range (M,M + dM), which are formed at the redshif t interval 
(z c , z c + dz c) and observed at redshift z (Sasaki 1994 ; Chiu & 
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Ostriker 2000) 

N(M, z, z c )AMAz c = N M (zc) 
D(z c ) 



D(zc) 
D(z)a(M)J D{z c ) 

dZc dM 



D(z) H(z c )(l + z c ) 



(1) 
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where the overdot represents derivative with respect to time. The 
parameter Nm (z c )dM is the number of collapsed objects per unit 
comoving volume within a mass ran ge (M, M + dM) at re dshift 
z c , known as the PS mass function (Press & Schechter 1974), and 
S c is the critical overdensity for collapse, usually taken to be equal 
to 1 .69 for a matter dominated flat universe (f2 m = 1 ) . This param- 
eter is quite insensitive to cosmology a nd hence the same value ca n 
be used for all cosmological models (Eke, Cole, & Frenk 1996). 
The other parameters are: H(z) is the Hubble parameter, D(z) is 
the growth factor for linear perturbations and a(M) is the rms mass 
fluctuation at a mass scale M. 

At early epochs, (i.e., redshifts larger than the reionization red- 
shift, z > Zre), the lower mass cutoff for the haloes which can 
host star formation will be decided by the cooling efficiency of the 
baryons. Indeed, the absence of heavier elements and lower free- 
zout molecular fraction of hydrogen (/h 2 — 10~ 6 ) in the primor- 
dial gas favours pure atomic cooling of the gas. This process can 
only make the gas cool upto 10 4 K. However, if one can some- 
how increase H2 (and probably HD) content of the gas, the final 
temperature can be reduced which can lead to the formation of 
cold gas condensations within the low mass haloes (T egmark et al. 
1997). Density enhancement during the collapse of the haloes can 
increase the molecular fraction upto 10~ 4 in central parts of the 
haloes (Haiman, Abel, & Rees 200C). However, the Lyman and 
Werner band photons that are pro duced by luminous objects can 



easily destroy these H2 molecules. Haiman, Abel, & Rees (200C ) 
have shown that the haloes smaller than T V1I — 10 ' K cannot cool 
even in the absence of any Lyman and Werner band flux. In what 
follows we consider models with T viT = 300, 10 3 and 10 4 K sepa- 
rately, which essentially covers the whole range of possibilities. 

After the universe has been reionized (z < z rc ), the lower 
mass cutoff for the haloes which are able to host star formation is 
set by the photoionization temperature (T ~ 10 4 K). It is known 
from simulations that the photoionizing background suppresses 
galaxy formation within haloes with circular velocities less than 35 
km s _1 , while the mass of cooled baryons is reduced by 50% for 
haloe^jvMi^^ircu^ar^eJadtifi^ & Weinberg 

1996). However, the exact value will depend on the intensity and 
spectral shape of the ionizing background radiation. We consider 
this as a free parameter spanning a range of 35 km s" 1 to 75 km s~ 
for the circular velocity cutoff (i>c Ut ). Note that Vc Ut is taken to be 
time-independent - hence, it does not corresp ond to the line ar the- 
ory Jeans mass (which decreases with time) ( |Gnedin 200C ). Also, 
in our model, the change from the neutral to reionized phase hap- 
pens abruptly at z = z rc . However, in realistic models one needs 
to consider the transition region where considerable fraction of gas 
remains in the ionized phase. This will lead to an extra suppression 
in the low mass haloes even in the pre-overlapping era. 

We assume that the SFR, of a halo of mass M which is formed 



at redshift z c , is given by ( Eggen^ Lynden-Bell, & Sandage 1962 



Cen & Ostriker 1 992| ; piiedin 1996| ) 



M S f(M,z,z c 



t 2 



x exp 



t(z) - t(z c ) 



tdyn 



(2) 



where esf is a efficiency parameter for the conversion of gas into 
stars. The function t(z) gives the age of the universe at redshift 2; 
thus, t(z) — t(z c ) is the age of the collapsed halo at z. We have 
assumed here that the duration of star formation activity in a given 
halo extends over a dynamical time-scale tdyn- The exponential de- 



crease of the SFR comes from the assumption that it is proportional 
to the mass of the cold gas (= total baryonic mass minus the mass 
already gone into stars) inside the halo (Schmidt's law). We use 
data obtained from rest-UV luminosities of galaxies, compiled an d 
corrected for extinction by Somerville, Primack, & Faber (2001), 



to normalize esf. Physically, this quantity depends on the metallic- 
ity of the gas and the effect of feedback due to outflows and local 
radiation field. For simplicity, we assume this to be independent of 
redshift for most of this article. Now we can write the cosmic SFR 
per unit comoving volume at a redshift z, 

p SF (z) = H dz c r dM'M SF (M',z,z c )xN(M',z,z c ).(3) 

Jz J M 

The lower mass cutoff M at a given epoch is decided by the Jeans 
criteria as explained above. 

Assuming the formation GRBs is closely related to the for- 
mation of stars, we get the birthrate of GRBs at redshift z per unit 
comoving volume as 



/Psf(z) 



(4) 



where / is a efficiency factor which links the formation of stars to 
that of GRBs. It is the number of GRBs per unit mass of forming 
stars, hence it depends on the fraction of mass contained in (high 
mass) stars which are potential progenitors of the GRBs. This im- 
plies that / is crucially dependent on the slope of the stellar initial 
mass function (IMF). 

Given the formation rate, it is straightforward to obtain the rate 
of GRB events per unit redshift range expected over the whole sky 
at present: 

^_«,,» ™ 

dz 1 + z dz 

where the factor (1 + 2) is due to the time dilation between z and 
the present epoch, and dV(z) is the comoving volume element 



dV(z) = 4.7TC 



H{z){l 



-dz 



(6) 



The parameter dh{z) is the luminosity distance. Note that the dif- 
ference between <j>{z) and diV GRB /dz is purely geometrical in na- 
ture. The number of GRBs per unit redshift range observed over a 
time At bs is then given by 



diV dD, . _ . dN^ H 
■3— = — At ohs $(z) — — 
dz 4.7T dz 



(7) 



where dSl/47r is the mean beaming factor and $(2) is the weight 
factor due to the limited sensitivity of the detector, because of 
which, only brightest bursts will be observed at higher redshifts. 



3 RESULTS AND DISCUSSION: 

The analysis presented here is based on LCDM model with the pa- 
rameters tt m = 0.3, Ov = 0.7, h = 0.65, n b h 2 = 0.02, as = 
0.93, n — 1. We consider three epochs of reionization (z rc ), three 
virial temperatures (T v i r ) upto which the gas can cool before reion- 
ization (i.e z > Zre) and three circular velocity cutoff (v^ ut ) for 
epoch after reionization. For a given set of model parameters, 
changing the reionization epoch is equivalent to changing the es- 
cape fraction of Lyman limit photons from the haloes. The SFR 
density obtained from our model is plotted as a function of red- 
shift in Fig |IJ In the panel (a) we plot the results for three epochs 
of reionization for T v i r = 10 4 K(solid curves) and 300 K (dashed 
curves) considering v£ ut — 35 km s _1 . 
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Figure 1. The SFR density as a function of redshift. In the panel (a) we plot the results for three epochs of reionization (z rc = 6, 9 and 12) for T v ; r = 10 4 
K(solid curves) and 300 K (dashed curves) considering u£ ut = 35 km s — 1 . In panel (b) the results are plotted for v^ ut = 35, 50 and 75 km s _1 considering 
T v i T = 10 4 K and z IC = 6. The curves are normalized using the extinction corrected data points from Somerville, Primack, & Faber (2001). 



Our models consistently reproduce the observed points in- 
ferred from rest-UV luminosities of galaxies at z < 5 for esF of 
0.1. Indeed the maxima in the curve occurs in the pre-reionization 
era. This is mainly due to the substantial contribution from the low 
mass haloes where the star formation is not suppressed by t he pho- 
toionization (Barkana & Loeb 2000, Bromm & Loeb 2002). How- 



ever, as expected, the SFR density depends strongly on the epoch 
of reionization and the nature of the cooling (hence T v j r ). Inter- 
esting point to note is that across the reionization epoch there is 
upto 3 fold increase in the SFR density. This should be considered 
as a strict upper limit as is achieved when the Lyman and Werner 
band photons are completely absent (hence T v i r = 300 K). Addi- 
tional difference of factor two could be introduced across z le by 
increasing the t£ ut to 75 km s - , thereby suppressing star forma- 
tion in more low-mass haloes in the post-recombination era, which 
is shown in panel (b). 

3.1 Dark GRBs: 

The optical after glows are only de tected in one third of the well 
local ized GRBs (pynbo et al. 2001 ; Lazzati, Covino, & Ghisellini 
2002; |Reichart & Yost 200l| ). The GRBs without afterglows are 
called dark GRBs. The optical limits obtained in the case of dark 
afterglows are fainter than the detected ones. This suggests that 
the dark GRBs can not be a ccounted for by the failu re to image 
deeply enough quickly (see Reichart & Price (20021 and refer- 



ences therein). The reason for these missing optical afterglows is 
attributed either (i) to the extinction in the host galaxy, or (ii) to the 
fact that some of the GRBs lie beyond the reionization epoch and 
the neutral IGM absorbs the afterglow. Indeed high- z QSO obser- 
vations show complete Gunn-Peterson a bsorption for z > 6 QSOs 
(|Fan et al. 200 1|; bjorgovski et al. 200lh. 



Our model predicts the redshift distribution of GRBs for a 
given /. Realistically, time evolution of the mass function of stars 
(for example triggered by the evolving metallicity) and the change 
in mean properties of the interstellar medium (ISM) will introduce 
time dependency for /. However, here we assume / to be constant. 



To predict the observed redshift distribution, one needs to take 
into account the fact that only the brightest bursts will be observed 
at higher redshifts due to limited detector sensitivity. We take this 
into account by assuming a luminosity function of the form 



L 



exp 



Lo 
' L 



(8) 



and then calculating the fraction of GRBs observed, <&(z), at a par- 
ticular z using the relation 



JT m , nW ^)^ 

f~1>(L)dL 



(9) 



Here, L is the "isotropic equivalent" intrinsic burst luminosity in 
the e nergy band 30 — 2000 keV (as defined by Porciani & Madau 
(2001)) and L m i n (z) is the corresponding minimum intrinsic lu- 
minosity that can be observed with the detector. We consider two 
separate cases - Case A: ip(L), and hence, Lo is independent of z, 
and Case B: ip(L) evolves with redshift, with the evolution given 
by Lo ~ (1 + z) 1 _^as_suggesterJ^y_Lloyd-Ronning, Fryer, & 
Ramirez-Ruiz (2002). 

For a given detector sensitivity, we can determine L m i n (z) 
from the relation 



Lmin(z) = P 



2000 kcV 
30 keV 



ES(E)dE 



(1 



x r(l+z)E n 
■V J(l + 2 )S m 



S{E)AE 



(10) 



where P is the minimum observed photon flux (sensitivity) at the 
detector in the energy band [_E m i n , Bmax] and S(E) is the rest 
frame GRB energy spectrum. S(E) is taken to be a broken pow- 
erlaw with a low-energy spectral index a, a high- energy spectral 
index f3 and a break energy Et (Band et al. 1992). In this work, 
we take a = —1 and = — 2.25, which are the mean values 
measured b y Greece ct al. (2000|), and Eh _ 511 keV (P orciani & 
Madau 2001). 

The minimum intrinsic luminosity at z that can be observed by 
the detector depends on z in two ways: (i) the decrease in the flux 
due to distance and (ii) K-correction. For the luminosity function, 
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Figure 2. Powerlaw spectrum of GRB afterglow at different redshifts (la- 
beled at emission wavelength of the Lyman-a line) modified by Gunn- 
Peterson effect of neutral IGM. The effect of damping wings is not con- 
sidered here. We also plot the window function of the R, I and J bands. This 
figure illustrates the fact that afterglow will be invisible in R and I band for 
z > 6.8 and 8.8 respectively. 



we take 7 _ —2.5, Lq 



3.2 x 10 51 ergs s 1 (which are obtained 



by 3 orciani & Madau (2001 ) using the log N — log P relation for 
BATSE bursts) and [E mia ,E max ] = [50,300] keV. We consider 
values of P corresponding to two detectors, namely, P — 0.2 
photons c m' 2 s" 1 (BATSE) and P = 0.04 photons cm" 2 s" 1 



(Swift) ( |Lamb & Reichart 2000| ). 

A similar analysis needs to be carried out for the afterglows 
too. According to the simplest afterglow model of W ijers, Rees, 
& Meszaros (1997), the spectrum follows a power-law form when 
observed in the optical and X-ray wave-bands about 1 day after 
the burst. A simple analysis shows that the afterglow of the faintest 
detectable GRB will always be brighter than 23rd magnitude within 
1 day after the burst in both I and R band. Hence, if a substantial 
fraction of GRB afterglows are not observed in I and R bands within 
1 day after the burst, it must be because of radiative transfer effects 
and not because of cosmological effects. 

We can now predict the percentage fraction of detectable 
GRBs at redshifts greater than 2, called F GRB (z). From Fig |J 
it is clear that due to high opacity of the IGM neither the afterglows 
nor the corresponding host galaxies will be detected in the R and 
I-band if the GRBs are located at z > 6.8 and 8.8 respectively. The 
fraction of R and I band dropout afterglows for a range of model 
parameters and detector sensitivities are given in the fourth to the 
seventh columns of the Table |I] (F GHB (R) and F GRB (1)). The 
results clearly indicate that even in the optimistic case (with lumi- 
nosity evolution and high sensitivity of the detectors), only 38% 
(25%)of the afterglows will not be detected in the R (I) band be- 
cause of the extinction due to the IGM opacity. It is interesting to 
note that even if the Hi optical depth of the ISM in the host galaxy 
is as high as 10 23 cm" 2 , due to the redshift effect, X-ray absorption 
produced by this gas will not affect our visibility of this source in 
the soft X-ray band (0.5 to 2 keV) whenever it is detected in the 
hard X-rays. However such a gas will produce damped Lyman-a 
absorption. 



Table 1. The predicted range in the percentage fraction of dark GRBs due 
to Gunn-Peterson effect, taking into account two different models for the 
luminosity evolution (Case A and B). The lower and upper limits of F are 
obtained with detector sensitivities similar to that of BATSE and Swift 



respectively. 




..cut 

u c 


T ■ 

± vir 


F GRB (R) 


F GHB (I) 




km s" 1 


K 


Case A 


CaseB 


Case A 


CaseB 


6 


75 


10 4 


0.8-3.9 


15.5-25.6 


0.2-0.9 


5.6-11.4 


6 


50 


10 4 


0.8-3.5 


14.1-23.7 


0.2-0.8 


5.1-10.5 


6 


35 


10 4 


0.7-3.3 


13.3-22.5 


0.1-0.7 


4.8-10.0 


6 


35 


10 3 


1.1-5.0 


20.0-32.9 


0.3-1.4 


8.8-17.7 


6 


35 


300 


1.3-5.7 


22.7-36.9 


0.4-1.7 


10.7-21.1 


9 


35 


10 4 


0.6-2.7 


12.3-21.8 


0.1-0.7 


5.2-11.1 


9 


35 


10 3 


0.9-4.1 


19.1-33.3 


0.3-1.5 


10.1-20.7 


9 


35 


300 


1.0-4.7 


22.1-37.9 


0.4-1.8 


12.4-24.9 


12 


35 


10 4 


0.4-1.9 


8.7-16.1 


0.0-0.5 


3.7-8.4 


12 


35 


10 3 


0.5-2.3 


12.3-24.0 


0.2-0.9 


7.3-16.7 


12 


35 


300 


0.5-2.5 


14.3-28.0 


0.2-1.1 


9.2-20.9 



3.2 Redshift distribution of GRBs: 

We next compare the differential redshift distribution of GRBs with 
that obtai ned using the luminosity - variability (L — V) c orrelation 
of GRBs ( Lloyd-Ronning, Fryer, & Ramirez-Ruiz 2002]; S chaefer, 
Deng, & Band 2001). The variability V of a GRB gives a mea- 
sure of the rms fluctuation of the burst light curve around the mean 
value. For the GRBs with known redshifts, it was seen that there 
is a correlation between the luminosity and the variability of the 

(Fenimore & 

+0.48 

0.54 



light curve, of the form L oc V", with a = 3.35_ 



Ramirez-Ruiz 2000), 3.3+n o (freichart et al. 200 1| ) and 1.57 



( Lloyd-Ronning, Fryer, & Ramirez-Ruiz 2002| ). The difference be- 
tween various results is because of the precise definition of variabil- 
ity used. Since the number of GRBs with known redshifts is small 
(~ 20), the scatter in the (L — V) correlation is very large. 

From this observed correlation, one can, in principle, obtain 
the the redshift dis tribution of GRBs. Lloyd-Ronning, Fryer, & 
Ramirez-Ruiz (2002) use a sample of 220 bursts to obtain the in- 
trinsic cumulative redshift distribution of GRBs JV b s (< 1 + z )- 
Their analysis, that corrects for the truncation effects in the sam- 
ple due to limited detector sensitivities, gives the intrinsic redshift 
distribution of the GRBs. We distributed this data into logarithmic 
redshift bins and calculated the slope [which is essentially the dif- 
ferential distribution dN b B /dz\ for each bin. The error in each bin 
is mainly contributed by the scatter in the L — V correlation. 

The quantity predicted by our model dN/dz (given by eq (Q), 
but without the factor <&(z) in the right hand side) contains a nor- 
malization factor /Aiobsdf2/47r, which is fixed by fitting to the 
data using \ 2 -minimization technique. In Fig^, we show the com- 
parison between our predictions and th e data obtained by L loyd- 
Ronning, Fryer, & Ramirez-Ruiz (2002) using the L — V corre- 
lation. Most of our models are broadly consistent with the data 
mainly because of large error bars. It is interesting to note that mod- 
els with T v i r < 10 3 K deviate a lot from the mean observed points; 
but large errors, mainly from the large scatter in the L — V re- 
lation, do not exclude these models. If the errors in the observed 
distribution reduce, one will be in a position to constrain the epoch 
of reionization and the contribution of molecular cooling (and other 
feedback mechanisms) in the formation of first generation of galax- 
ies. 
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Figure 3. The redshift distribution of GRBs for different model parameters. 
The data points with error bars are from the observed L — V correlation 
(Lloyd-Ronning, Fryer, and Ramirez-Ruiz 2002). The error bars along the 
y-axis represent the 1<t error due to the scatter in the observed L — V cor- 
relation, while the error bars along the x-axis show the bin size. The results 
from our model are shown for z rc = 6 and 9, and for three different values 
of T vir - namely, 10 4 K (solid curves), 10 3 K (dotted curves) and 300 K 
(dashed curves). 



3.3 Discussions and Conclusions 

In this article we have modelled the redshift distribution of GRBs 
assuming that they follow the cosmic star formation history. The 
main conclusions are summarized below: 

(i) A peak in the SFR before reionization is seen because of sub- 
stantial contribution from the low mass haloes where the star for- 
mation is not suppress ed b y the photoionization ( also noticed by 
Barkana & Loeb (200C|) and|Bromm & Loeb (200^)). We show that 



the redshift distribution of the SFR density depends on the reion- 
ization epoch and the minimum virial temperature (T v i r ) of star 
forming haloes. The minimum virial temperature is set by the na- 
ture of the cooling - if H2 cooling is effective, haloes with much 
lower virial temperatures can participate in the star formation, 
(ii) Our results clearly indicate that even in the optimistic case (with 
luminosity evolution and high sensitivity of the detectors), only 
38% (25%)of the afterglows will not be detected in the R (I) band 
because of the extinction due to the IGM opacity. This means that a 
substantial fraction of optically dark GR Bs (> 66%) originate be- 
cause of effects such as dust extinction ( Waxman & Prairie 2000| 



Ramirez-Ruiz, Trentham, & Blain 2002; Berger et al. 2002) 



(iii) Although our predictions vary a lot as we change our model 
parameters, we cannot constrain the parameter space effectively 
because of large scatter in observations. One needs an improved 
observed redshift distribution to make real progress in understand- 
ing the physical processes in epochs prior to the reionization. 

(iv) In general, the efficiency parameters esf and / depend of red- 
shift. The star forming efficiency esF evolves because of its de- 
pendence on the metallicity of the collapsing gas and the effect of 
feedback due to outflows and local radiation field. Similarly, the pa- 
rameter / ev olves as the s tellar IMF becomes less top-heavy with 
decreasing z ( Larson 199E ). With improved observational data, one 
can use our model to constrain the evolution of these parameters. 

One would naively expect that the presence of damping wing 
of Hi absorption could be useful for probing the neutral IGM. How- 



ever, in a simple picture of GRBs following supemovae, the pres- 
ence of Hll region and damped Lyman-a absorption because of the 
host galaxy will complicate the matter. Thus, we believe that the 
redshift distribution with less scatter may be a better way to probe 
the reionization. 

We thank N. M. Lloyd-Ronning for providing us with data 
on the cumulative redshift distribution of GRBs. We also thank T. 
Padmanabhan for useful comments. T.R.C. is supported by the Uni- 
versity Grants Commission, India. 
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